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Abstract 

A parameter-free calculation for J/ty suppression in pA collisions, based on 
the dipole description, is confronted with the new data from the PHENIX 
experiment. Achieving good agreement, we employed this model predicting 
the contribution of initial state interactions (ISI) to J/^f suppression in A A 
collisions. Such a transition from pA to AA is not straightforward, since 
involves specific effects of double color filtering and boosting of the saturation 
scale. Relying on this refined ISI contribution, we updated the previous 
analysis of RHIC data on J/^> production in Cu-Cu and Au-Au collisions 
at yfs = 200 GeV, and determined the transport coefficient of the created 
dense medium at go = 0.6 GeV 2 / fm. Nuclear effects for J /\I/ production 
at the LHC are predicted using the transport coefficient go = 0.8 GeV 2 / fm, 
extracted from data on suppression of high-p^ hadrons in central lead-lead 
collisions at y/s = 2.76 TeV. Our analysis covers only direct J/\l/ production, 
while data may also include the feed-down from decay of heavier states and 
i?-mesons. 
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1. Introduction 

Charmonium production in heavy ion collisions originates from two sources, 
which can be labeled as initial (ISI) and final (FSI) state interactions. The 
former are also called sometimes cold nuclear matter effects, and are related 
to the creation of a cc pair and its subsequent propagation through the fast 
moving (in the collision c. m.) matter, consisting of nucleons and their rem- 
nants. The second stage, FSI, is related to the propagation and attenuation 
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of the J/^f through the (probably) dense matter created in the collision. This 
stage is usually the main point of interest, since it is considered as a probe 
for the properties of the created medium. 

The two stages of collision are characterized by very different time scales, 
so no interference between them is possible. It has been always an important 
and challenging task to separate the effects originated from different stages 
of interaction. Essential information about the ISI stage can be withdrawn 
from data on J/ty production in pA collisions, where it cannot be mixed up 
with the effects of FSI. Although a model-independent extrapolation of ISI 
from pA to AA collisions is impossible [lj, this is a good playground for the 
models. 

In Section [2] we perform parameter free calculations for J/^> production 
in pA collisions at \/s = 200 GeV and different rapidities. The results agree 
well with the recently released data from the PHENIX experiment at RHIC. 
Predictions for the LHC energies are provided as well. 

In Section [3] we extend the model to heavy ion collisions, including the 
specific effects of double color filtering and boosting of the saturation scales. 
These effects do not exist in pA collisions. The pr-dependence of the nuclear 
ratio Raa is predicted at the energies of RHIC and LHC. This first stage of 
calculations cannot yet be compared to any data, since the effects of FSI are 
to be added. 

in Section @] we calculate the additional suppression factor related to at- 
tenuation of the produced J/^f due to FSI with the created medium. The 
key observations here are the short formation time of the J/\l/ wave function, 
and the relation between the break-up cross section of a cc dipole and the 
transport coefficient of the medium. If the latter is not known, the calcula- 
tions cannot be done in a parameter-free way, but can be compared with data 
attempting to determine the transport coefficient. We found the parameter 
q = 0.6 GeV 2 / fm for gold-gold collisions at y/s = 200 GeV. For the energies 
of LHC we rely on the analysis of the data on high-p^ hadron suppression in 
lead-lead collisions at y/s = 2.76 TeV, which resulted in q Q = 0.8 GeV 2 / fm 
and provided predictions for J/ty suppression. 

Our results are displayed as function of pr, as well as pr-integrated. 

2. J/\l/ production in pA collisions 

Any time scale is subject to Lorentz transformation and at high energies 
it may be long. In particular, the time scale for charm quark production, 
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which is very short, t* ~ l/m c in the proper reference frame, may become 
longer than the nuclear radius in the nuclear rest frame, 
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This condition is well satisfied at y/s — 200 GeV, where l c = 12 fm at y = 0, 
and rises as e y at forward rapidities. Therefore, a cc pair is not produced 
momentarily inside the nucleus, but it propagates through the whole nucleus. 
So, inclusion of the break up of the produced cc dipole, as is usually done, 
is not correct. Notice that this should not be a debatable issue, if one wants 
to comply with Lorentz invariance. 

The long life time of a cc in the projectile hadron significantly enhances 
nuclear suppression. One can see that on the example of a very heavy nucleus. 
If J/ty were suppressed only by inelastic processes, leading to break up of 
the cc dipole produced on mass shell, the A-dependence of the cross section 
would be A 2 / 3 . If however, the cc pair propagated through the whole nucleus, 
the A-dependence would be A 1 / 3 . 

The break up cross section for a cc dipole is usually assumed to be un- 
known and is fitted to the same data which is to be explained. However, 
this cross section is well known from precise measurements at HERA. The 
proton structure function F 2 (x,Q 2 ), measured in a wide range of Bjorken x 
and virtuality Q 2 , provides a sensitive probe for the dipole (break up) cross 
section, dependent on the dipole size and energy. Many parametrizations 
of the dipole break up cross section are currently available in the literature. 
One of the most popular results from a global analysis of HERA data gives 
a break up cross section parametrized in the saturated shape 0]. In what 
follows we rely on this result. 

In the regime of t c 3> Ra, when the dipole size is "frozen" during propa- 
gation through the nucleus, the nuclear suppression factor was calculated in 



1|, |3J as, 
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RpA = J d 2 b J dz \S pA (b,z)\' 
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where 



S pA (b,z) 



d 2 r T Wc C {r T ) exp 



--(Tccg(r T )T_(b, z) - ^cr 5c (r T )T + (b,z) 



(3) 
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Here W 5c (r T ) oc K {m c r T ) r\ ^j/^(r T ) is the distribution function for the 
dipole size r^; <Jc Cg (r T ) = ^a Sc (r T /2) — |cr gc (r T ) is a three-body (gcc) dipole 
cross section; T_(6, z) = f_ dz' pA{b, z')\ T + (b,z) = T A (b) —T-(b,z), and 
T A {b) =T_(6,oo). 

We employ the phenomenological dipole cross section o" ec (rT,x) fitted to 
HERA data in Q. Here 

x= vmrw) e ., t (4) 

v s 

where y is the rapidity, and the mean invariant cc mass squared was fixed at 
the value corresponding to the color singlet model, (M| c ) = 2Mj^. 

Here we improve the accuracy of the calculations of [l|, relying on a 
more realistic evaluation of the J/\& mean radius, (rfyj/y = § (fj/^), where 

J(r 2 j^) = 0.42 fm 0, Q. In Fig. Q] we compare the results with the new 

data for J/^ suppression in deuteron-gold collisions at = 200 GeV, 
and provide predictions for LHC. This calculation includes gluon shadowing, 
based on the global NLO analysis of DIS data Q (see details in although 
the correction is quite small. The agreement with data is rather good, taking 
into account that this is not a usual fit to the data to be explained, but is the 
first parameter free calculation. Notice that the extension of our calculations 
down to y < —1 should be taken with precaution, since the conditions for 
using the frozen dipole approximation, as well as dipole phenomenology itself, 
are not fulfilled. 

We also calculated the nuclear modification of the pr-dependence, known 
as Cronin effect. Following jsl, 0] we introduced into the integral over b in 
Eq. (J2J) a pr-dependent factor, normalized to unity, 

R p a(pt, b) = }f T } * pA (h . (l + 7l + w 2 . g . m ) " 6 ■ (5) 

(p 2 T ) + A pA (b) \ 6(p 2 T )J V 6[(Pt) + A pA(b)} J 

Here (pj) is the mean transverse momentum squared for J/^> produced i n pp 



collisions. For the energy dependence we use the parametrization from [10 
(Pt) = [—2.4 + 0.6 Ins] GeV 2 . A. pA (b) = (pt)pA~ (Pt)pp is nuclear broadening 



for a J/f produced at impact parameter b. It was calculated in [Hj in good 
agreement with the available data for J broadening, 

A pA (b) = ^C q (E J/ < D )T A (b), (6) 
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Figure 1: Data [6] for the p^-integrated nuclear suppression factor RdAu{y) f° r J/^ pro- 
duced in d-Au collisions with rapidity y at i/s = 200 GeV. The upper solid curve presents 
the result of Eqs. ©-(13]). It includes the small correction for gluon shadowing The 
lower solid (dashed) curve shows predictions for LHC at ^Js — 5.5 TeV, including (exclud- 
ing) gluon shadowing. 



where the factor C q (Ej/y) was calculated in [U 



The results of the calculation of the nuclear ratio R p a{pt) are depicted 
in Fig. |2] as function of px, at rapidities y = 0, 2 and at <Js = 200 GeV and 
5.5 TeV. 

Notice that these predictions are not reliable above px > 5 GeV, where 
data for pp collisions are lacking, and the px dependent factor Eq. flS} cannot 
be justified. An alternative shape of the px dependence, suggested in (sf, 
leads to similar results at px % 5GeV, but a smaller Cronin effect at larger 
p T . 



3. ISI effects in nuclear collisions 

With this model for the nuclear effects related to ISI, tested in pA colli- 
sions, one can move on to heavy ion collisions, and predict their ISI effects, 
which is desperately needed in order to single out from data the FSI sup- 
pression caused by the dense medium produced in the collision. 
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Figure 2: Prediction for the p^-distributions of Jj^l produced with rapidities y = 0,2 in 
p-Au collisions at ^fs = 200 GeV and in p-Pb collisions at ^/s = 5.5 TeV. 



As was pointed out in [1| the transition from pA to AA collisions is not 
straightforward and is model dependent. There are several effects specific for 
nucleus-nucleus collisions, and one of them is double color filtering. Namely, 
color filtering of cc dipoles propagating through the first nucleus makes the 
mean dipole separation smaller, so the second nucleus becomes more trans- 
parent compared to attenuation of a cc dipole produced in pA collisions. The 
ISI suppression factor in a collision of nuclei A and B, including double color 
filtering was derived in [I], 



RAB(b,r) 



;i+A A + A+)(l + A+ + A B ) 
(A+ - A- A )(A + B - A B ) "* L(l+A+ + A+)(l + A^ + A£) 



In 
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where b and f are the impact parameters of the collision and of the produced 
J/ty respectively. The other notations are, 
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2rl(x A ) 
32r 2 (x A ) 
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A « = Mb)*™-^' <10) 

As was mentioned above, we rely on the parametrization of the dipole cross 
section 0], with parameters o"o = 23.03 mb and r (x) = 0.4 x (x/xg) ' 144 fm, 
and x = 3.04 x 10~ 4 . The value of xa is given by Eq. (jlj), and xb = %a e 2y ■ 
The mean separation squared (rf.) of the produced cc dipole was evaluated 
in Q at (r|) = 0.045 fm 2 . 

The ISI stage also modifies the px distribution of produced J/^s, as was 
discussed in the previous section. The full effect of ISI, including attenu- 
ation and the modified px dependence, has the form of the ratio Eq. ([7j) 
multiplied by the normalized factor, Eq. flSJ), with the replacement A p A(b) =>■ 
&-ab(P, f) = A p a(t) + A P B(b — f). The pr-dependent nuclear ratio for J/ty 
produced with y = in central Cu-Cu and Au-Au collisions at \fs = 
200 GeV; and in Pb-Pb at = 2.76 and 5.5 TeV, is depicted in Fig. H 
These calculations include all the ISI effects, but exclude FSI. 

We also calculated the effects of boosting the saturation scale in nuclear 



collisions [12J], which leads to an increase of broadening A^a(&, t). The en- 
hancement factor Kaa(P,t) was calculated for J/ty production in pj. We 
plotted the results by dashed curves, since this effect is still under debate 
and has not been confirmed by data so far. 

4. FSI suppression and the transport coefficient 

The nuclear effects related to ISI shown in Fig. H] present the results of 
parameter-free calculations, which are, however, model dependent. Although 
the model was tested in Section [2] comparing with d-A data, the transition 
from pA to AA collisions is very nontrivial [l| and involves more modeling, as 
was discussed in the previous section. Nevertheless, relying on these results 
for the ISI contribution, one can try to pin down the effects of FSI, which 
probe the properties of the medium created in nuclear collisions. Such an 
analysis of RHIC data for J/^f production in heavy ion collisions at RHIC 
was performed in 0,0]. Here we update the analysis relying on a more refined 
calculation of the ISI effects, and also make predictions for the LHC. 

The key point of the analysis [sj for the FSI stage is the relation between 



the dipole cross section and broadening 13j, [l4|, which allows to express J/ty 
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Figure 3: Nuclear ratio for J/^i> produced in central (b = 0) AA collisions including only 
the ISI and excluding FSI effects. Solid lines correspond to (from the top to bottom) Cu- 
Cu and Au-Au at y/s = 200 GeV, Pb-Pb at y/s = 2.76 and 5.5 TeV. The dashed curves 
are different from the solid ones by inclusion of the effect of boosted saturation scale (see 
text for explanations). 

attenuation in terms of the transport coefficient of the medium, 



An important observation also employed here, is the shortness of the forma- 
tion time of the J/\I/ wave function. This is why the transverse cc dipole 



The transport coefficient q is defined as the magnitude of broadening ex- 
perienced by a quark propagating through a path length 1 fm in the medium 
[1 51 ] . With the usual assumptions that the initial medium density at t — to is 
proportional to the number of participants n part , and the density is diluting 
with time as 1/t, we arrive at 




(12) 






(13) 



t n. 



''part 
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where the parameter q is the maximal value of q, for the medium produced 
at t = to in central collision at b = r = 0. As previously, b and f are the 
impact parameters of collision and position of the quark, which is propagating 
through the medium in the transverse direction (in the collision cm.). 

Combining the J/^f attenuation factor caused by FSLEq. ffl2l) . with 
the ISI effects depicted in Fig. [31 as is described in jsl, 0], we arrive at 
the final nuclear modification of the J/ty production rate, which can be 
now compared with data, as is shown in Fig. HI The parameter go = 
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Figure 4: RHIC data 3, 17, IH for the nuclear modification factor Raa for J/\P produc- 
tion in central collisions of Cu-Cu (closed circles and squares) and Au-Au (open circles) at 
y = and yfs = 200 GeV. The curves are the same as in Fig. |3l but corrected for the FSI 
effects, which are calculated with the transport coefficient for gold-gold qo = 0.6 GeV 2 / fm 
adjusted to the data. Raa at the LHC energies is predicted with go = 0.8 GeV 2 / fm, 
extracted in [l9| from data for nuclear quenching of high-pr hadrons 



0.6 GeV 2 / fm for gold-cold collisions was adjusted to RHIC data lfj, [l7|, Il8 
shown by solid and open circles for Cu-Cu and Au-Au respectively. This 
parameter was rescaled for copper-copper according to Eq. (Tl3l) as (go)cu = 
Kart(0,0)/< rt (0,0)](g )A u = 0.38 GeV 2 / fm. 

In order to predict Raa at the energies of LHC, one needs to know the cor- 
responding values of %. This parameter was found to be % = 0.8 GeV 2 / fm 
in the analysis 19j of the recently published data 20] on suppression of high- 
Pt hadrons in central lead- lead collisions at y — and y/s = 2.76 TeV. We 
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use this value of g at both energies y/s = 2.76 and 5.5 TeV, since do not 
expect a significant variation of go in this energy range. The results for the 
LHC are depicted in Fig. H] by the bottom pair of solid curves. 

Dotted curves show the effect of mutual boosting of the saturated scales in 
colliding nuclei. We see that at y/s = 200 GeV this correction is quite small 
and is visible only at > 5 GeV, where calculations, as was mentioned 
above, are not reliable anyway. At the energies of the LHC the correction is 
larger and may be detected if data were sufficiently precise. 

We also calculated the nuclear modification factor for the J/^f production 
cross section integrated over px at y = and different energies. The results 
are presented in Table [TJ 

Table 1: The pT-integrated nuclear ratio for J/ty production at y — in central collisions 
of different nuclei at the energies of RHIC and LHC. 



Cu-Cu 




= 200 GeV 


Raa 


= 0.54 


Au-Au 




= 200 GeV 


Raa 


= 0.34 


Pb-Pb 


V~s 


= 2.76 TeV 


Raa 


= 0.18 


Pb-Pb 




= 5.50 TeV 


Raa 


= 0.16 



5. Summary 

Summarizing, we tested parameter- free calculations for J/ty suppression 
in pA collisions, based on the dfpole description 0, [lj], with the new data 
from the PHENIX experiment [6J]. Encouraged by the good agreement, we 
employed this model predicting the ISI contribution to J/^f production in 
AA collisions. Such a transition from pA to AA involves novel effects of dou- 
ble color filtering and boosting the saturation scale [H, which were included. 
The previous analysis [g| of RHIC data on J '/\& production in Cu-Cu and 
Au-Au collisions at y/s = 200 GeV is updated, and the transport coefficient 
of the created dense medium is found to be go — 0.6 GeV 2 / fm. The esti- 
mated accuracy of this analysis is about 30%. To predict nuclear effects for 
J/ty production in Pb-Pb at the LHC we rely on the value of the trans- 
port coefficient go — 0-8 GeV 2 / fm, determined from the analysis of data on 
suppression of high-pr hadron production in central lead-lead collisions at 



v^ = 2.76 TeV [20 . 
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Finally, a word of caution is in order: the calculations are done for direct 
production of J/^f, while currently available data may be significantly con- 
taminated with a feed-down from decays of x, an d -B-mesons. This may 
significantly disturb the magnitude of the observed nuclear effects, because 
heavy charmonia apparently are more suppressed by FSI than J/^f, while 
S-mesons are expected to be much less affected by the ISI and FSI effects 
21 . 
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